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Abstract 

Recent ineleastic electron-proton scattering experiments have led to rather accurate 
values for the N — > A transition quadrupole moment Qn->A- The experimental results 
imply a prolate (cigar-shaped) intrinsic deformation of the nucleon. The nonsphericity 
of the proton's charge distribution might be seen in the spectrum of atomic hydrogen. 
The possibilities and limitations for determining the geometric shape of the nucleon in an 
atomic physics experiment are discussed. 

1 Introduction 

Elastic electron-proton scattering experiments and atomic spectroscopy measurements have 
shown that the positive charge of the proton is distributed over a finite volume, and that the 
charge radius of the proton is about r p = 0.9 fm [1]. Furthermore, the proton's radial charge 
density p(r) could be determined from the Fourier transform of the charge form factor measured 
in elastic electron-proton scattering. However, the proton charge distribution most certainly 
has an angular dependence, i.e., p(r) = p(r, 6). This angular dependence of the proton charge 
density is not directly accessible in elastic scattering. 

Inelastic photon and electron scattering experiments provide further details about the inner 
structure of the nucleon. In particular they show that the proton has a spectrum of excited 
states. When photons of sufficient energy impinge on a hydrogen target, the proton maybe 
excited to its lowest lying excited state with spin 3/2, called the A + (1232), where the number 
in parentheses is the energy (mass) of this state in units of MeV. The A + resonance decays 
quickly via the strong interaction into the proton ground state of 939 MeV and the lightest 
strongly interacting particle, the pion, with a mass of about 140 MeV. The lifetime of the A + 
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is of the order of 10~ 23 s, i.e., the time it takes for a gluon to travel across the proton radius of 
1Q-15 m ^nis corr esponds to a line width of about 100 MeV. 

From a multipole analysis of the electromagnetic p — > A + transition amplitude we have 
infered that the proton charge density deviates from spherical symmetry. In fact, recent exper- 
imental determinations of the p — > A + transition quadrupole moment [2, 3, 4] provide evidence 
for a prolate (cigar-shaped) intrinsic deformation of the proton [5]. 

2 Structure of the nucleon 
2.1 Quark model 

The structure of the nucleon 2 ground state and its excitation spectrum is explained in the quark 
model [6]. In the quark model the proton and neutron consist of three spin 1/2 quarks which 
come in two charge states, called up and down quarks, i.e., p{uud) and n(ddu). Here, u and d 
denote the up-quark with electric charge e u = 2/3 and down-quark with charge e<2 = —1/3. The 
first excited states of the proton and neutron, the A + (uud) and A°(ddu) with spin J = 3/2, 
have the same quark content as the p and n with spin J = 1/2. Thus we expect the properties 
of the A to be closely related to those of the N ground state. Formally, N and A properties are 
related because there exists a higher symmetry than isospin, called SU(6) spin-flavor symmetry, 
which combines the N isospin doublet with J = 1/2, and the A isospin quartet with J = 3/2 
together with other states in the same 56-dimensional spin-flavor supermultiplet. 

The strong interaction between the quarks is modelled as a long range harmonic oscillator 
confinement force, which guarantees that the colored quarks cannot be observed as free particles. 
In order to explain the spectrum in greater detail, for example, the mass splitting between 
the iV(939) with spin 1/2 and the A(1232) with spin 3/2 (see Fig. 1), the long-range spin- 
independent confinement forces must be supplemented by short-range spin-dependent terms. 

Spin-dependent forces arise, e.g., from one-gluon exchange between quarks. The one-gluon 

2 The proton and neutron, which have nearly the same mass of about 939 MeV can (with respect to strong 
interactions) be regarded as one particle, called nucleon 7V(939), that exists in two charge states. Similarly, 
the first excited state of the nucleon, the A(1232) exists in four charge states A - , A , A + , A ++ , which are 
energetically almost degenerate, i.e., all four states have a mass of about 1232 MeV. Mathematically, this 
symmetry is decribed by the SU(2) isospin group, whose generators are, in analogy to the spin, the Pauli 
matrices. The N has isospin 1/2 and the A has isospin 3/2. The quarks building up the N and A have 
approximately the same mass and thus have isospin 1/2. 
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Figure 1: The hyperfine splitting between the iV(939) and A(1232) due to the spin-dependent 
interaction between quarks. This splitting is the nucleon physics analogue of the 21 cm line 
between the hyperfine states F = and F — 1 in atomic hydrogen. In nucleon spectroscopy, 
an excited state is usually denoted by the name of the resonance followed by its energy [MeV] 
in parentheses. 

exchange potential (OGEP) is analogous to the Fermi-Breit interaction between the electron 
and proton in atomic hydrogen, except for a factor A« • Xj, which reflects that quarks carry color 
charges, and the replacement of the photon-electron coupling constant a by the strong gluon- 
quark coupling a s . The most important terms of the one-gluon exchange potential (spin-orbit 
terms are omitted) are 

y° GEP (^ ',-)=t* • A ^ i 1 + h • 5(r) - i^p (3<T * • ^ • • *4 (1) 

Here, r = i-j — Tj is the distance between the i-th and j-th quark; r = r/r is the corresponding 
unit vector, Oi is the spin operator (Pauli matrices), and Xi is the color operator (Gell-Mann 
matrices). The constituent quark mass of the up and down quarks is denoted as m q . Usually 
one chooses m q = Mn/3 where Mn is the nucleon mass. In the entire paper we use natural 
units H — c — 1. 

The 1/r term is called color-Coulomb interaction, and the ai ■ Oj term is refered to as 
color-magnetic interaction. The matrix element of the operator Xi ■ Xj between color singlet 
states, such as the iV and A, is readily evaluated and gives a factor -8/3. Therefore, the color 
magnetic interaction is repulsive between quark pairs coupled to spin 1 and attractive between 
quark pairs coupled to spin 0. 

The mass difference between the N and A masses of 293 MeV is almost completely deter- 
mined by the color- magnetic interaction in Eq.(l); the N — A mass difference is often refered 
to as hyperfine splitting (see Fig. 1) in analogy to the 10" 6 eV splitting between the F = and 
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F — 1 hyperfine states in the hydrogen atom ground state. The huge difference of fourteen or- 
ders of magnitude between the hyperfine splitting in the nucleon and hydrogen atom is mainly 
due to the different size of both systems 3 . 

2.2 The electromagnetic N — > A transition 

If the A resonance is produced in an electromagnetic process (see Fig. 2), parity invariance 
and angular momentum conservation restrict the iV — > A excitation to magnetic dipole (Ml), 
electric quadrupole (E2), and charge or (Coulomb) quadrupole (C2) transitions. Although 
the E2 and C2 amplitudes are small (about 1/40 of the dominant magnetic dipole transition 
amplitude) they are nevertheless important because their nonzeroness provides evidence for a 
deviation of the nucleon charge distribution from spherical symmetry [2, 3, 4, 5]. In the quark 
model, the A(1232) is obtained from the iV(939) ground state either by flipping the spin of a 
single quark (Ml transition) or by flipping the spins of two quarks (E2 or C2 transition) [7]. 
In the case of the quadrupole transitions there is also a small contribution coming from e.g., 
lifting a single quark from an S state in the nucleon into an orbitally excited D state in the A. 
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Figure 2: The excitation of the A(1232) resonance by a virtual photon 7 of momentum Q 
is described by the three electromagnetic transition form factors G^* A (Q 2 ), G^ A (Q 2 ), and 
Gq^ A (Q 2 )- Their contribution to the inelastic electron-nucleon cross section can be obtained 
by analyzing the angular distribution of the decay pions in coincidence with the scattered 
electron. In the limit Q — > 0, the cross section is described by the magnetic dipole moment 
[J> N—*A and the charge quadrupole moment Qn^a- 

With the quark model it is also possible to calculate how much the proton charge dis- 
tribution deviates from spherical symmetry and which degrees of freedom are responsible for 

3 The matrix element of the <5-function in Eq.(l) is proportional to 1/R 3 , where R is a measure of the spatial 
extension of the system, i.e., the quark core radius b w 0.5 fm in the case of the nucleon, and the Bohr radius 
ao ~ 0.5 A in the case of the hydrogen atom. 
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this deformation [5]. A multipole expansion of the nucleon charge density operator p up to 
quadrupole terms leads to the following invariants in spin-isospin space 



3 3 

i i<j 



2a i ■ Oj — (3cr i Z o j z — (Ji ■ a j) 



(2) 



where <jj Z is the z-component of the Pauli spin matrix of quark i, and e { = |(1 +r iz ) is the 
quark charge where Tj 2 is the third component of the Pauli isospin matrix. The constants 
A and B in front of the one- and two-quark terms parametrize the orbital and color matrix 
elements so that p is only an operator in spin-isospin space. It is important to note that the one- 
body operator (the A term) in Eq.(2) represents the valence quark degrees of freedom, whereas 
the two-body operator (the B term) even though it also acts on quark variables provides an 
effective description of the quark-antiquark degrees of freedom in the physical N and A [5, 8, 7]. 
Three-quark currents are omitted here, but were taken into account in Ref. [9]. 

After evaluating Eq.(2) between quark model spin-isospin wave functions the following re- 
lations between the p, n, and A + charge radii, and between the p — > A + transition quadrupole 
moment Q p ^a+ and the neutron charge radius r 2 were derived [7] 

r l - r l+ = r l, Qp^a+ = ^= r 2 n , (3) 

where r 2 p and r\+ are the proton and A + charge radii. Inserting the experimental neutron charge 
radius [10] r 2 n = —0.113(3) fm 2 in the second relation of Eq.(3) we obtain Q p ^a+ — —0.08 fm 2 . 
This agrees well with recent determinations of Q p ^a+ from inelastic electron-proton scattering 
data which yield Q p ^ A +(exp) = -0.0846(33) fm 2 [3], and Q p ^ A +(exp) = -0.108(9) fm 2 [4]. 
The second relation is the zero momentum transfer limit of a more general relation between the 
iV — > A charge quadrupole transition form factor Gq^* a (Q 2 ) and the elastic neutron charge 
form factor G^(Q 2 ) 



G^ A (Q 2 ) = - 3 -^GUQ% (4) 



which is experimentally satisfied for a wide range of momentum transfers [11]. 

These relations are a consequence of the underlying SU(6) spin-flavor symmetry and its 
breaking by the spin-dependent two-body terms in Eq.(2). As we have noted before, the spin- 
spin interaction in Eq.(l) is repulsive between quark pairs in a spin 1 state and attractive in 



quark pairs with spin 0. This explains why the A + , which contains only quark pairs coupled 
to spin 1 is heavier than the proton. Similarly the electromagnetic counterparts of these spin- 
dependent terms in Eq.(2) explain why r\+ > r^, why the neutron charge radius is negative, 
and why the neutron has a prolate intrinsic deformation. This can be qualitatively understood 
as follows. The two down quarks in the neutron are always in a spin 1 state 4 . Consequently, 
the spin-spin force pushes them further apart than an up-down quark pair. This results in an 
elongated (prolate) charge distribution with the up quark in the middle, and at the same time 
in a negative neutron charge radius. On the other hand, in the neutral A, where all quark 
pairs are in a spin 1 state, there is an equal distance between up-down and down-down quark 
pairs. This corresponds to an equilateral triangle (oblate) configuration of the charges and a 
zero charge radius of the neutral A (see Fig. 3). 



Figure 3: Qualitative picture of the neutron (left) and A (right) charge distributions in the 
quark model. Although the deformation is depicted here as residing in the valence quark 
distribution, in reality a spherical quark core is surrounded by a deformed cloud of quarks and 
antiquarks [5]. A similar picture is obtained for the p and A + by interchanging u and d quarks. 



4 In order to satisfy the Pauli principle the total three-quark wave function, which is the direct product of 
orbital, spin-isospin, and color wave functions, must be completely antisymmetric under permutation of any 
two quarks. Because the color singlet N and A wave functions are antisymmetric in color space, the product 
of orbital and spin-isospin wave functions must be symmetric. For orbitally symmetric ground state S-wave 
functions this in turn means that the spin-isospin wave functions must be completely symmetric. The latter 
can only be achieved if the two down quarks in the neutron, which are necessarily in an isospin 1 state, are 
simultaneously in a spin 1 state. 
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3 Shape of the nucleon 

3.1 Intrinsic quadrupole moment of the nucleon 

In order to learn something about the shape of a spatially extended particle one has to determine 
its intrinsic quadrupole moment. The intrinsic quadrupole moment of a nucleus 

Q = Jd'rp(v)(3z 2 -r 2 ), (5) 

is defined with respect to the body-fixed frame. If the charge density is concentrated along 
the z-direction (symmetry axis of the particle), the term proportional to 3z 2 dominates, Qo is 
positive, and the particle is prolate (cigar-shaped). If the charge density is concentrated in the 
equatorial plane perpendicular to z, the term proportional to r 2 prevails, Q is negative, and 
the particle is oblate (pancake-shaped). 

The intrinsic quadrupole moment Qo must be distinguished from the spectroscopic 
quadrupole moment Q measured in the laboratory frame. In the collective nuclear model 
[12], the relation between the observable spectroscopic quadrupole moment Q and the intrinsic 
quadrupole moment Qo is 

3K 2 -/(/+l) 
Q ~ (J+l)(2/ + 3) Qo > (6) 

where / is the total spin of the nucleus, and K is the projection of / onto the z'-axis in the 

body fixed frame (symmetry axis of the nucleus) 5 . The factor between Q and Qo is due to the 

quantum mechanical precession of the deformed charge distribution with body-fixed symmetry 

axis z' around the lab frame z-axis along I 2 . Therefore, in the lab frame one does not measure 

the intrinsic quadrupole moment directly but only its projection onto the lab frame quantization 

axis. As can be seen from Eq.(6) a spin 1 = nucleus does not have a spectroscopic quadrupole 

moment Q even if its intrinsic quadrupole moment Qo is different from zero. This is also the 

case for the spin 1=1/2 proton. 

Recently, we have determined the intrinsic quadrupole moment of the proton and A + in the 

quark model [5] and found 

Ql_= ~rl Q A+ = rl (7) 

5 In the following the nuclear spin is denoted as I and the total angular momentum of the electrons is denoted 
by J. 



7 



i.e., the intrinsic quadrupole moment of the proton is given by the negative of the neutron 
charge radius and therefore positive, whereas the intrinsic quadrupole moment of the A + is 
negative. This corresponds to a prolate proton and an oblate A + deformation, consistent with 
the qualitative explanation given above. 

The same result, namely a connection between the neutron charge radius r\ and the intrinsic 
quadrupole moment of the proton Q p is also obtained in the pion cloud model (see Fig. 4). In 
the pion cloud model, the nucleon consists of a spherically symmetric bare nucleon (quark core) 
surrounded by a pion with orbital angular momentum 1 = 1. For example, the neutron can be 
viewed as being composed of a bare proton surrounded by a negative pion. The nonspherical 
pion wave function leads to prolate intrinsic deformation of the nucleon. For further details see 
Ref. [5]. 




Figure 4: Intrinsic quadrupole deformation of the nucleon (left) and A (right) in the pion cloud 
model. In the N the p-wave pion cloud is concentrated along the polar (symmetry) axis, with 
maximum probability of finding the pion at the poles. This leads to a prolate deformation. 
In the A, the pion cloud is concentrated in the equatorial plane producing an oblate intrinsic 
deformation (from Ref. [5]). 

Using for the nucleon the simple model of a homogeneously charged rotational ellipsoid, we 
can give Qo a geometric interpretation in terms of the half-axes of the ellipsoid. In classical 
electrodynamics the quadrupole moment of a rotational ellipsoid with charge Z, major axis a 
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along, and minor axis b perpendicular to the symmetry axis is given by 

Q = ^(a 2 -b 2 ) = ^ZR 2 5, (8) 

with the deformation parameter 5 = 2(a — b)/(a + b) and the mean radius R = (a + b)/2. 
Inserting the quark model result Qq = — r 2 = 0.113 fm 2 on the left-hand side of Eq.(8), and for 



R the equivalent radius R = ^5/3r p , where r p is the proton charge radius, we obtain a nucleon 
deformation parameter 5^ = 0.11. This corresponds to a ratio of major to minor semi-axes 
a/b = 1.11. For the deformation parameter of the A we find 5a = —0.09 and a half-axis ratio 
a/b = 0.91. 

3.2 Intrinsic quadrupole charge density of the nucleon 

To proceed, we decompose the proton charge form factor in two terms, a term resulting from a 
spherically symmetric charge distribution, and a second term due to the intrinsic quadrupole 
deformation of the actual charge density 

G P c (Q 2 ) = G p vol (Q 2 )- 1 -Q 2 G p def (Q"). (9) 

For the spherically symmetric part we take the usual dipole form factor G D (Q 2 ). Concerning 
the intrinsic quadrupole charge density, we employ the relation between the N — * A charge 
quadrupole transition form factor and the elastic neutron charge form factor in Eq.(4), and the 
relation between the intrinsic quadrupole moment of the nucleon and the neutron charge radius 
Eq.(7). We then obtain the following expressions 

G p vol (Q 2 ) = G D (Q 2 ) = (1 + Q 2 /AT 2 , G p def (Q 2 ) = -V2G^ A (Q 2 ) = ^ 2 GUQ 2 ) 
G P c(Q 2 ) = G p vol (Q 2 )-G n c (Q 2 ), (10) 

with Gj e j(0) = Qq. Thus, the deviation of the neutron charge form factor from zero and the 
deviation of the nucleon's charge distribution from spherical symmetry are closely related. 

According to Eq.(9) the proton charge density in coordinate space is decomposed in two 
parts, a spherically symmetric volume term p p Q i(r), and a nonspherical charge density denoted 

as Pdef(r) 

P p {r) = p P vol (r) + lv 2 p p def (r)=pl ol (r)-p n (r), (11) 
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where the last equality follows from Eq.(lO). These coordinate space expressions show that the 
deviation of the proton charge density from spherical symmetry is described by the neutron 
charge density p n (r). To see the effect of the nucleon's nonsphericity explicity, we use a two 
parameter fit [13] for the elastic neutron charge form factor 

(IT 

G n c (Q 2 ) = -^ TT ^G D (Q 2 ), (12) 

where r = Q 2 /(4M^) and A 2 = 0.71 GeV 2 . We then obtain for the Fourier transform of 
Eq.(12) 

1 A 4 rr) 2 / p-mr -kr -kr A 

™ - V> ( m> — - A2 — - C" 2 - A2 »< — - < 13 > 

where m =: 2M^/\fd. The neutron structure parameters a and d have been determined from 
the lowest moments of the experimental neutron charge form factor. They are given by the 
neutron charge radius r 2 , and the fourth moment r\ [14]. Similarly, the Fourier transform of 
the intrinsic charge quadrupole form factor G p de j(Q 2 ) yields 

p , , 2 A 4 m 2 fe~ mr e- Ar m 2 - A 2 _ Ar \ 

with Jd 3 r fJ% e f(r) = Qo = —r 2 . Finally, the spherically symmetric part of the proton charge 
density is obtained from the Fourier transform of the dipole form factor 

Ploi(r) = ^e^. (15) 

The decompostion of the proton charge density in a spherical volume part and a quadrupole 
deformation part is shown in Fig.5. 

4 Atomic spectroscopy and nucleon structure 

Before discussing the possibility of measuring the intrinsic quadrupole moment of the nucleon 
in an atomic physics experiment, I would like to recall that atomic spectroscopy provided 
the first clear evidence for nuclear quadrupole moments before these were seen with nuclear 
spectroscopic methods and nuclear scattering experiments. The quadrupole moment of the 
hydrogen isotope D was discovered using the molecular beam magnetic resonance method [15] . 
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Figure 5: Spherically symmetric charge density of the proton p p vo i{r) [fm~ 3 ] (broken curve); 
intrinsic quadrupole charge density of the nucleon p^ e f(r) [fm _1 ] (broken-dotted curve); neutron 
charge density p n (r) = — | V 2 Pde/( r ) (thin full curve); and the sum p p (r) = p v vol — p n (r) (thick 
full curve). As a result of the intrinsic quadrupole deformation of the nucleon positive charge 
is moved to the exterior region of the proton. This will lead to a larger proton charge radius 
and will be discussed in sect. 4.3. 

4.1 Discovery of nuclear quadrupole moments in atomic hyperfine 
spectra 

Early in 1935, when measuring the magnetic hyperfine structure of certain atomic spectral lines 
of europium isotopes 151 Eu and 153 Eu with nuclear spin 7 = 5/2 Schiiler and Schmidt found 
small but systematic deviations from the Lande interval rule [16]. The latter describes the 
magnetic hyperfine splitting of spectral lines due to the magnetic interaction energy between a 
nucleus with magnetic moment fij and spin I, and the magnetic moment of the atomic electrons 
of total angular momentum J 

AW hfs = A < I • J >= X -AC, A — ^H(0)/(IJ). (16) 

Here, A is a constant proportional to the magnetic field 77(0) produced by the atomic electrons 
at the site of the nucleus. Furthermore, C = F(F + 1) — 7(7 + 1) — 7(7 + 1) is the Casimir 
factor where F = I + J is the sum of the nuclear and electronic angular momenta. 

Schiiler and Schmidt noted that their experimental results could be explained if a deviation 
of the nuclear charge distribution from spherical symmetry is assumed. The correct quantum 
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mechanical interpretation was given by Casimir [17] who showed that the observed hyperfine 
splitting pattern could be reproduced by adding a second term on the right hand side of Eq.(16) 
of the form 

AW 1 AC I 1 q 3C(C+1)-4/(/ + 1)J(J+1) B _ n 7) 

where the quadrupole constant B is given by the product of the nuclear quadrupole moment Q 
and the electric field gradient due to the electrons at the site of the nucleus. The quadrupole 
moment in Eq.(17) is refered to as the spectroscopic quadrupole moment. From the measured 
hyperfine splitting AWhfs and the calculated electric field gradient the quadrupole moment 
Q could be determined. For the europium isotopes the result was Q( 151 Eu) = 150 fm 2 and 
Q( 153 Eu) = 320 fm 2 . 

For 7 = and 7 = 1/2, Eq.(6) predicts Q = and a vanishing quadrupole contribution 
to the hyperfine splitting AWhf s even if the intrinsic quadrupole moment Q ^ 0. Because of 
these selection rules it seems at first sight hopeless to obtain information on the shape of an 
7 = 1/2 nucleus, such as the nucleon, from an atomic physics experiment. However, there are 
other observables that are sensitive to the geometric shape of the nucleus. 

4.2 Isotope shifts and intrinsic quadrupole moments 

The frequency of a certain spectral line of a given element is slightly different for the various 
isotopes of this element. These isotope shifts arise even in the absence of nuclear electromagnetic 
moments because different isotopes have different mass and size. Heavier isotopes with bigger 
radii experience less Coulomb binding than lighter isotopes with smaller radii, because the 
attractive inner atomic Coulomb potential is cut off near the surface of nucleus. Consequently 
the spectral lines of the heavier isotope are slightly shifted to lower frequencies compared to 
the same lines of the lighter isotope. For s electrons, the isotope shift is given by the following 
expression 

5E IS = —\^ e (0)\ 2 5<r 2 >+5E mass + 5E pol , 
5 <r 2 > = 5 < r 2 > vol +5 < r 2 > def (18) 

where 5 < r 2 > is the change in the nuclear charge radius of the two isotopes considered. 
Here, 5 < r 2 > vo i and 5 < r 2 >def denote the volume and deformation contribution to the total 
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charge radius change. The electronic wave function at r = is denoted as ^/ e (0) = ^[Z/(na )] 3 , 
and a is the Bohr radius. The terms 5E mass and 5E po i are the mass and nuclear polarization 
contributions to the isotope shift. In the following we discuss only the size (volume) and shape 
(deformation) contributions to the isotope shift. 

If the radius change in Eq.(18) were a pure volume effect, we would expect the radius 
variation to be based on R = RqA 1 ^, where A is the number of nucleons and R = 1.2 
fm. This relation between nuclear mass and equivalent radius R = \J\ < r2 > follows from 
the nuclear liquid droplet model, which considers a nucleus as a spherical droplet of constant 
density. On the other hand, the deformation contribution to the charge radius change is related 
to the difference of intrinsic quadrupole moments Q of the considered isotopes. Thus, there 
can be a substantial radius increase if the charge distribution of one isotope is more deformed 
than the other even if their volumes are nearly the same. 

Already in 1934 Schiiler and Schmidt [18] found an anomalously large isotope shift between 
the spectral lines of the spin 1 = nuclei 150 Sm and 152 Sm, which was nearly twice as large 
as between neighboring isostope pairs, of mass number 148 — 150 and 152 — 154. This large 
shift could not be explained by a pure volume effect 6 . The explanation offered by Brix and 
Kopfermann in 1947 [19] is based on the following idea: (i) the anomalous isotope shift is 
connected with the jump of the quadrupole moments between 151 Eu and 153 Eu, and (ii) the 
1 = nuclei 150 Sm and 152 Sm have approximately the same intrinsic deformation as the 151 Eu 
and 153 Eu nuclei with I = 5/2. In other words, they assumed that the spectroscopic quadrupole 
moments observed in the europium isotopes were already present in the samarium isotopes. 
This seemed reasonable because the large quadrupole moments of the europium isotopes could 
not possibly be generated by the addition of a single proton to the corresponding samarium 
nuclei. With this assumption they could explain the observed large isotope shift between the 

two samarium isotopes. 

pure volume effect would give a radius change of 5 < r 2 >= 0.21 fm 2 , whereas the observed value was 
approximately S < r 2 >= 0.48 fm 2 . 
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Figure 6: Decomposition of the N and A charge radii into a volume and a deformation con- 
tribution. The intrinsic quadrupole deformation of the nucleon and A leads to an increase of 
the proton (left) and A + (right) charge radii. The broken lines correspond to the spherical 
volume contribution. The deformation contribution to the charge radii indicated by the thin 
ellipsoids is given by the negative of the neutron charge radius according to Eq.(19). The full 
circles represent the measured charge radii which contain the effect of averaging the deformed 
charge distribution over all directions in space. 

4.3 Estimate of the nucleon deformation contribution to hydrogen 
hyperfine structure 

Applying these ideas to the nucleon, we decompose the proton and A + charge radii in two 
terms, a spherically symmetric volume contribution and a deformation contribution due to the 
intrinsic quadrupole charge density (see Fig. 6). This corresponds to the decomposition in 
Eq.(ll). The deformation contribution to the proton and A + charge radii, which makes the 
charge radius bigger, is determined by the neutron charge radius 

2 2,2 2 2 , 2 

ry ry __l_ ry ry ry _|_ ry 

'p ~ ' p,vol ^ ' p,defi A+ A+,vol ' A+,def 

r 2 p = r 2 p + r 2 n + (-r 2 n ), r 2 A+ = r 2 p + (-r 2 n ). (19) 

If written in this way we see that the volume contribution to the proton charge radius is 
given by the isoscalar charge radius r 2 ol = rj s . The latter is defined in terms of the proton 
and neutron charge radii as r 2 IS =: r 2 + r 2 . The deformation contribution to the proton charge 
radius is given by the negative of the neutron charge radius r pde f = —r 2 . Similarly, the volume 
contribution to the A + charge radius is given by r\+ vol = r 2 , and the deformation contribution 
is (as in the case of the proton) given by r\ + de j = —r 2 . 

In the quark model the isoscalar charge radius is mainly determined by the quark core radius 
b 0.6 fm and the charge radius of the constituent quark r q ^ 0.6 fm 

r 2 IS = b 2 + r 2 . (20) 
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The quark core of the nucleoli and the constituent quarks themselves are assumed to be spher- 
ical. Concerning the quark core, this is supported by explict calculation of the P and D wave 
probabilities in the nucleon wave function, which are very small [8]. The deformation of the 
nucleon arises predominantly from the collective quark- ant iquark degrees of freedom described 
by the two-body operators of section 2.2. 

To estimate the energy shift of s electrons in atomic hydrogen which results from the 
deformation of the nucleon's charge density we use 5 < r 2 >= —r 2 in Eq.(18). For n = 1 
and a = 0.529 A we then obtain 

5E IS (dei) = 157 kHz. (21) 

Compared to other corrections in the hydrogen hyperfine structure [20] this seems to be a big 
effect. We emphasize that this deformation contribution is already included in the experimen- 
tal proton charge radius [21]. One possibility to isolate the deformation contribution to the 
proton charge radius is to measure hydrogen isotope shifts. The increasing accuracy of such 
experiments [22] will perhaps reveal some discrepancies between theory and experiment that 
arise from the intrinsic deformation of the hydrogen isotopes, and which manifest themselves 
as deformation contributions to their charge radii. In this connection we mention that the 
proton charge radius as determined in hydrogen Lamb shift measurements may be principally 
different from the one determined in elastic electron-proton scattering. Atomic physics exper- 
iments involve small electronic line widths and thus measure the electromagnetic interaction 
of the bound electron with the proton over long time scales. Thus they see a time average of 
the proton's deformed charge distribution. Because of this averaging effect they yield a bigger 
proton charge radius. On the other hand, high energy electrons involve comparatively short 
interaction times, which corresponds to taking a snapshot of the deformed charge distribution 
having a particular orientation in space and thus see a slightly smaller charge radius. 

5 Summary 

In summary, the electromagnetic excitation of the nucleon to its first excited state, the A 
resonance, has provided clear evidence that the charge distribution of the nucleon deviates 
from spherical symmetry. The intrinsic quadrupole moment of the nucleon, which is defined 
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with respect to a body-fixed coordinate system that co-rotates with the nucleoli, is a measure 
of the nucleon's quadrupole deformation. We have calculated the intrinsic quadrupole moment 
of the nucleon in the quark model and found that it is given by the neutron charge radius, 
implying a prolate deformation of the nucleon's charge distribution. More generally, we have 
suggested that the neutron charge density p n (r) is a measure of the intrinsic quadrupole charge 
density of the nucleon p p de f{r). 

Isotope shifts of atomic spectral lines have provided information on the intrinsic deforma- 
tion of spin 1/2 and spin nuclei, which do not have spectroscopic quadrupole moments due 
to angular momentum selection rules. The hydrogen spectrum can be measured with very high 
accuracy. It does not seem completely unlikely that some future experimental technique, per- 
haps involving muonic hydrogen, or hydrogen molecules, such as H 2 , HD, and D 2 , will reveal 
further nucleon structure details such as the spatial shape of the proton. 
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